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Abrikosov vortex lattice dynamics in a superconductor with weak defects is studied taking into 
account gyroscopic (Hall) properties. It is demonstrated that interaction of the moving lattice with 
weak defects results in appearance of the additional drag force, F, which is F(V) tx W at small 
velocities, while at high velocities F(V) tx Thus, the total drag force can be a nonmonotonic 

function of the lattice translational velocity. It leads to the velocity dependence of the Hall angle 
and a nonlinear current-voltage characteristic of the superconductor. Due to the Hall component in 
the lattice motion, the value [dF(V)/dV]~ 1 does not coincide with the differential mobility dV/dF e , 
where F e is the external force acting on the lattice. The estimates suggest that condition dV/dF e < 
can be met with difficulty in contrast to dV/dF < 0. The instability of the lattice motion regarding 
nonuniform perturbations has been found when the value of (-dF/dV) is greater than a certain 
combination of elastic modulii of the vortex lattice, while the more strict requirement dV/dF e < 
is still not fulfilled. 

PACS numbers: 74.20.De, 74.25. Fy, 74.60.-w 



I. INTRODUCTION 

Vortex dynamics in superconductors is still a contro- 
versial matter in physics of superconductivity. As a rule 
dynamics of Abrikospv vortices was considered in the 
overdamped regime,™ when the viscosity is strong and 
it namely determines lattice dynamics. The Hall force 
was generally neglected in conventional superconductors, 
where viscosity is high. Interest in the Hall force has 
quickened in the past few years with fabrication of high. 
T c superclean single crystals with the large Hall angle.H 
Moreover, among other things, Hall properties are man- 
ifested in special collective modes existing in the vortex 
system. These modes were observed in resonant experi- 
ments with 4 HcJ3 The interest to these modes as applied 
to high-T c superconductors has increased considerably ip. 
recent years-in connection with experiments in BSCCOu 
and YBCOcil compounds in which resonant effects were 
observed. In YBCO system resonances occur most likely 
because of cyclotronic vortex modes, while the resonant 
effects in BSCCO are associated either with Josephson 
plasma modes emerging them owing to Josephson cou- 
pling of adjacent CuC>2 layeraj or with the vortex modes 
excitations.Q 

The vortex modes may manifest themselves not only 
in resonances but also in dissipation of the moving vortex 
lattice caused by excitation of these modes. This effect is 
owing to irreversible transfer of the translational vortex 
motion energy into elementary bosonic excitations - the 
quanta of vortex modes in the presence of the vortex- 
defects interaction. For the single vortex it has been 
demonstrated that excitation of vortex modes results in 
an anomalous behavior of the friction force, F v , which is 
divergent at V — > OB The next step is study of vortex 



modes excited in the moving vortex lattice. 

In the present work we investigate vortex lattice dy- 
namics in the presence of weak defects. We are inter- 
ested in an additional effect of dissipation induced by the 
vortex lattice motion and interaction with weak defects, 
which are not strong enough to pin the vortex lattice, 
but this motion results in excitation of normal vibra- 
tional modes and irreversible transfer of the interaction 
energy to these modes (Landau damping). 

Due to the additional dissipation effect nonlincarity 
of the lattice motion has been revealed and it has been 
demonstrated that the nonlinear regime may be even 
dominant one in lattice dynamics being compared with 
the bare linear dissipation. When the lattice is the sub- 
ject of the Hall force we predict appearance of a lattice 
motion instability. The Hall angle is found to be lattice 
velocity dependent and the higher velocity the larger Hall 
angle as far as it is saturated. 

The paper is composed of five sections as follows. In 
Sec. II we describe a microscopical model of vortex lattice 
dynamics with the Hall and viscous contributions. In Sec. 
Ill we study the steady motion of the vortex lattice and 
discuss the stability of this motion. Discussion of the 
results of the model is presented in Sec. IV. Finnaly, the 
summary is given in Sec. V. 



II. MODEL 

Let us consider dynamics of the vortex lattice (formed 
by Abrikosov vortices parallel to the z— axis). Let vor- 
tices in equilibrium to be placed in the points of an ideal 
lattice, I. Dynamics of the vortex lattice can be described 
on the basis of an effective equation for the 2D vector 
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u = ui(z, t), lying in the xy— plane and describing the dis- 
placement of the I - th vortex in the lattice (see Fig.l).0 
The following Lagrangian determines nondissipative dy- 
namics of u: 



dz 



H 
~2 



du 



Ui, 



l,y 



at 



dui. x 
' dt 



-W{ui} - U imp {u}. 
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Here H is the Hall constant of the single vortex per 
unit length, W — W{u{\ is the energy of the ideal (de- 
fectless) deformed lattice, Ui mp {ui} defines the energy of 
interaction of the vortex with defects (its structure will 
be considered below) . To describe dynamics in the linear 
approximation VF{u/} should have the following struc- 
ture: 



w ^i} = ^E/{ k 
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where 5ui : i> = ui(z,t) — ui'(z,t), U lk (l) are components of 
the force tensor which are expressed through the second 
derivatives of a well-known potential of the vortex- vortex 
interaction with respect to equilibrium vortex positions,el 
k is the energy per unit length. We neglect in the above 
expression of the Lagrangian the inertial term, which 
comprises the vortex mass because, as was demonstrated 
irJ3, the vortex mass plays a substantial role in the dissi- 
pation of the moving vortex lattice only at high enough 
velocities. In the framework of a macroscopic approach, 
i.e. if ui varies insignificantly over a distance of the order 
of the lattice constant a v , the last expression (^) goes 
over to the standard expression as with the theory of 
elasticity, 
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where 011,044, Cee are the elastic modulii of the lattice, 
C44 can be expressed through k, C44 = n/a 2 . 

Each vortex moving in the medium experiences the ac- 
tion of the drag force and the rate of its energy dissipation 
is proportional to (dui/dt) 2 . For this reason, we choose 
the dissipative function in the form: 
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dui 
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dz, 



(4) 



where j3 is the dissipation coefficient per unit length of 
the vortex. When the motion of the vortex lattice is 



translational and steady (ui(z,t) = u/(0) + Vt, where V 
is the translational vortex lattice velocity) the dissipative 
function yields the friction force per unit length of the 
vortex, F = ~V(Q/V 2 ) = -V(3. 

In the case of "weak" defects the microscopical origin 
of defects is not essential.0 We determine interaction be- 
tween the vortex lattice and defects in the crystal as a 
variation of a local critical temperature, T c , depending 
on coordinates. This can be taken into account by in- 
troducing a coordinate-dependent coefficient [clq + F(r)] 
into the Ginzburg-Landau expansion!] We assume that 
the distribution of the order parameter in the vortex 
does not change during the vortex motion and is di 
scribed by a known function f(r±), where = x 2 +y 
Thus, regarding the displacement of the l-th vortex, we 

can write the vortex lattice energy ^ |^o| 2 /(nj_)) where 

1 

[(x -l x - ui x f + (y-l y - u iy ) 2 ] 1/2 and * is the 



order parameter. In this case, the energy £/j mp associated 
with crystal inhomogeneities can be written in the terms 
of the function f(r±) in the form of a functional of the 
vortex lattice displacement: 



E / dxdydz f(r±)F[ x + l x + ui x ,y + ly + ui y , z],(5) 



where x = x - l x - ui x , y = y - l y - ui y . 

Finally, one arrives at the equation for the displace- 
ment of the I— th vortex, ui(z). It is presented in the 
form: 



H z x 
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(6) 



Here the left hand side is the dynamical term governed 
by the Hall (gyroscopic) constant, the terms on the right 
hand side are forces acting on the I— th vortex because 
due to interaction with other vortices in the ideal lattice, 
friction and the presence of defects. 

Using the Eq. (|5]), the force caused by defects can be 
written in the form of the Fourier expansion in x and y: 



Ft 



l,imp 



^^/(gx^^e^'+^'e^'W , (7) 



ft ^ 



where ft = L x L y L z is the superconductor volume, J-(q) 
is the Fourier transform of the function F(r), which de- 
fines inhomogeneities in the system, see Eq.(^), q± = 
(q x , q y , 0), f(q±) is the vortex form-factor. 

In the absence of dissipation and defects the equation 
(^) can be simplified and because of translational invari- 
ance the force tensor components U tk depend on I — V 
only and do not depend on z. Then, the equation of 
motion can be presented as it is done in the lattice dy- 
namics theory, i.e. by means of— the Fourier transform 
of the force tensor components.E3 By using the Bloch 
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theorem it is possible to present the equation of mo- 
tion in the normal coordinates Uq >x (y) considering that 

tilMv) = S %x(») ' exp (iqj_l + iq z z\ . 

The equation of motion acquires the following form: 



Hii q 
-Hu„ 



Cyy(o) ' u q,y + ^xy\Q) ' u q.x'i 



(8) 



where Cs, 



the problem more complicated as equations are still lin- 
ear. However, the presence of defects results in nonlinear 
equations as the interaction energy (@) and force (^|) de- 
pend nonlinearly on u. A general solution of this nonlin- 
ear equation (^) cannot be found. In line withQ we can 
carry out a complete analysis assuming that deviations 
from the uniform translational and steady motion of vor- 
tices due to defects are small. Putting u = e x Vt + u(r, t) 
and linearizing (pi) in respect to u, we obtain the expres- 



J2 U lk (l) ■ exp(iql ) is Fourier representa- s i on for F irnp without ui(z): 



tion of the force tensor U lk , see (|). Here and below 
dot denote the time derivative. The Fourier transform 
of ui is discrete in the xy~ plane and continuous along 
the z— axis, the quantity q± has the sense of a quasimo- 
mentum, and q z is a regular linear momentum, (we put 
h = 1). The quasimomentum q± there is within the first 
Brillouin zone, \q\ < l/a„, whereas q z is only restricted 
by the physical reason, i.e. q z < l/£, where £ is the coher- 
ence length of the superconductor. Though these values 
q± and q z , have somewhat different physical meanings, 
below we will use a macroscopical approach and retain 
these symbols. 

The variables comprise the time dependent mul- 
tiplier, Hq~ tx exp[±iuj(q)i\, where uj(q) is periodic with 
regard to q±, tu(q + g) — uj(q), g is the vector of the re- 
ciprocal lattice. Then, one can derive the dispersion law 
using new parameters C\ — C xx — A, C2 = C yy +A, 2A = 



(C X X Cyy) + \j (CxX ^yy) + 



4C 2 • 



(9) 



At large |g| ~ l/a v the dispersion law of vortex lat- 
tice oscillations becomes periodical in the q± plane, with 
a period determined by the reciprocal lattice constant 
and the maximal value q Xj y = l/£ since the coefficients 
C X x(yy) an d C xy are periodical functions of q±. Near 
the center of the first Brillouin zone, |g| <C l/a„, the 
dispersion law ([)]) can be presented via elastic moduli, 
Ci = a 2 v {caql + c 66 q 2 ± ) ,C 2 = a 2 v (c i4 q z + Cuq^): 



H 



(cAAq 2 + ce6q±) (caq 2 + cng^) 



(10) 



which has been also obtained in Rcfs.^,|TT|. This law cor- 
responds to a low frequency gapless mode in the limit of 
small q z and q±. 

From here on we will employ a Debye-like model, with 
approximation of the dispersion law by the long wave 
equation (10) and with taking into account periodicity 
in the q Xl q v directions. The dispersion law for the vor- 
tex lattice is schematically drawn in Fig. 2. The Debye 
frequency for u>(q±), lud, is marked in Fig. 2 by a doted 
line. 

We have concerned lattice oscillations without dissi- 
pation and defects and found the dispersion law of nor- 
mal modes. The dissipation consideration does not make 
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(ii) 



Thus, the lattice motion through the nonuniform 
medium results in excitation of small oscillations of vor- 
tices in the lattice (vortex lattice normal modes), which 
are described by u. It was discussed by the authors of 
Ref. |l2] in the terms of "the lattice melting". Never- 
theless, it is essential for us that oscillations bring into 
appearence of the additional contribution to dissipation 
of the lattice, AQ = ((3/2) J2i J (diii/dt) 2 dz, and in turn 
into the additional drag force, F(V), 



p {v) =-v!m, F{v) =^m. (12) 



The general equation for u q (compare with (g)) reads 
1 



D 
1 

D 



[(Cyy + ifaxV) fx - (C Xy - lHq x V) fy] J (13) 
[(Cyy + l/3q x V) fy ~ (C X y + lHq x V) f X ] , 



where C lk = C ik (q), D = (C\ + i0q x V) ■ (C 2 + i(3q x V) 
(Hq x V) 2 , f x , v is related to @ 



fx, y = ^qx, y f(q±)f(q)e 



-iq^Vt 



(14) 



It is clearly seen that the condition D = at (3 = 
defines the dispersion law found in (||), if one would 
replace u> — > q x V . It is a typical evidence of resonant 
character of excitation of normal modes because of the 
lattice uniform motion. 

To analyze specifically a contribution of these small 
oscillations to dissipation of the moving vortex lattice we 
assume that inhomogeneity is caused by a system of point 
defects whose size is smaller than the radius of the vortex 
core. In this case, the function F(x,y,z) in (|J) can be 
written as the sum of Dirac delta functions: 



F(r) = } aS(f- r a ). 



(15) 



Here r a is the coordinate of the a— th defect and a 
characterizes the intensity of interaction of the defect. 
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The force acting on the I— th vortex with regard to the 
model (15) is written as (0) with substitution 



(16) 



By substituting (|T^) with the taking into account Jl6| ) 
in the dissipative function, employing the transformation 



dt 



: . _> / d 3 q(q x V) 2 u q u q , 



(17) 



where star denotes the complex conjugation, and carry- 
ing out averaging over defects with the help of the rela- 
tion ^exp(i<7r a ) = Ni m p5 qi o, where Ni mp is the number 

a 

of defects, 5 qi o is the Kronecker symbol, we can find the 
additional dissipation of the moving lattice. It may be 
associated with the drag force F(V). In the framework 
of the model (15) one can obtain the equation for the 
additional drag force F(V) per unit length of the vortex 



F(v)-r pVa2 

F(V)-C imp — 



d 3 q-q 2 J 2 (q±) 



q\ [q x V 2 (l3 2 + H 2 ) + C 2 ] + (Ci + C 2 )(q x VA - q^fGy - C 2 - A) 

[CiG 2 - (/? 2 + w){q x vy]> + p 2 (q x vy(c 1 + c 2 y 



(18) 



where Ci mp = Ni mp /Q is the defect concentration. 

Because of the factor (3V appearing in (O) in front 
of the integral one could assume a linear velocity depen- 
dence of the drag force, F(V) oc j3V, where (3 is some 
constant. Then, it seems that a consideration of vibra- 
tions of vortices associated with interaction with defects 
results only in a small correction to the bare friction force, 
Fbare = — @V . However, the integral in (|l^), as well as 
for the case of the single vortex,H contains singularities, 
therefore, the additional relaxation channel can become 
significant and even predominant. The most vivid exam- 
ple of this fact is that the additional dissipation can be 
finite in the limit of zero bare dissipation^ F(V) is finite 
at — ► 0, as F bare -> 0. 

At first glance, this appears as paradoxical. However, 
such a colissionless damping (Landau damping) stem- 
ming from the energy transfer from one mode to an- 
other appears in many branches of physics. For instance, 
flux- flow dynamics caused by excitations of low frequency 
fermionic modes has been proposed in Ref. [l^ and the 
finite friction force independent of the bare friction is 
predicted there. The friction force finite at V — ► and 
caused by excitation of bending oscillation of the domain 
wall in ferromagnets with microscopic defects was pre- 
dicted in Ref. As it will be demonstrated below here 
F(V) oc V s , S < 1 at V -> and F(V) > F bare at any (3 
and small enough velocities. 

This behavior 
The expression (18) 

(3/ {(3 2 + G 2 [q, (q x V) 2 ] }, where the function G is such 
that the condition G [q, (q x V) 2 ] = with the substitu- 
tion q x V — > oj defines the frequencies of normal modes of 
vortex vibrations in the lattice (|l0|). For (3 — ► the ex- 
pression ( |l8| ) is transformed into the <5-function, and be- 
comes ttS\G \q, (q x V) 2 ]}. After simple transformations 
this expression can be reduced to the 5-function of the 



can be explained as following, 
contains (3 in the combination 



type S [(q x V) 2 — u 2 {q)] , where ui(q) is the frequency of 
normal mode (|l0|). Since the equation q x V = u>(q) pos- 
sess a solution at any finite velocity the value of Q(V) 
can be finite even at (3 — > 0. 

In the case f3 — > 0, it is obvious that the additional 
dissipation becomes the main source of dissipation in the 
system. This may be realised if the friction force has a 
dependence like F oc V s , S < 1, as it was revealed in the 
case of the single vortex. 

There is an alternative explanation of the additional 
dissipation appearance. The modes excitation in the 
limit of small dissipation can be described on the ba- 
sis of momentum and energy conservation. Let us go 
over to a reference frame moving with vortices. In this 
reference frame defects move at a velocity —V — Ve x 
parallel to the x— axis and can transfer the momentum q 
to the vortex lattice only simultaneously with the energy 
<}V — qxV ■ This momentum is redistributed between 
the lattice as a whole and an elementary excitation, and 
the energy is transferred to elementary excitations only. 
In particular, in the case of the single vortex the vortex 
as a whole acquires momentum perpendicular to its axis 
(z— axis), while the wave propagating along the vortex 
obtains the z— axis momentum projection. 

Let us analyse the general expression ([18]). Note that 
integration has to be done not only inside the first Bril- 
louin zone. The integration is limited only by cut off q± 
over l/£, when the vortex form factor f(q±) 0. In the 
case of the vortex lattice, the lattice constant a v is usu- 
ally assumed to be larger than £. In order to resolve the 
vortex lattice problem one should go beyond the long- 
wave limit. Generally speaking, one has to integrate ( |l8| ) 
taking into account complicated periodical dependencies 
of Ci, C*2, A of the wave vectors, and a non-periodic 
dependence of f(q±), terms with q x V and so on. Nev- 
ertheless, the situation is simplified in limiting cases, for 
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example, at small or high velocities. 

First we concern the case of small dissipation, (3 -C H, 
for which the conservation law uj(q) — QxV defines an 
area that contributes to the integral in (|18|). To provide 
an estimation, we assume that en ~ cqq ~ c±, then 
luH = oJ;(c44(j£ + cj_g5.), and the condition u)(q) = q x V 
may be written as 



c 4 ACj_q 2 z + [c±q x - VH/2a v f 



i i 

c ±-% 



[VH/2a v ) 



(19) 



This expression indicates that at small V the size of the 
region in the qj_ plane, which contributes to the integral 
in ( p^| ) is proportional to V. If one considers that en ^ 
C66, then of course the above consideration remains true 
and only the expression (|l^) becomes more complicated. 
Thus, at V — > the region near q = gives rise to a 
very smalLcontribution to Q, which is proportional to 
V s , 6>lB 

However, if periodicity is taken into account, then the 
areas with u>(q) — > emerge near all points in the q± 
plane, which are equivalent to the origin. If one puts 
q± = g "+ q±, where qj_ is the wave vector within the first 
Brillouin zone and g is the wave vector of the reciprocal 
lattice, then these points are located at q± — > 0. It is 
worth to mention that values of q x V near these points 
are (g x + q x )V ~ 2nnV/a v . The corresponding area in 
the q± plane is now given by oj{q) = 2imV/a v . 

The value of g x is limited by the cutoff g x < l/£ or 
n < (a v /2T:^). Hence, the condition of the area smallness 
in the q^ plane with substitution of q x ~ g x ~ l/a v in 
q x V and with q^ in C\,Ci without changes takes the 
form: 



a 2 c 



(20) 



Thus, these regions can be small being compared to 
the size of the first Brillouin zone at velocity V -C V c> 
where V c — £c/H ~ where c is some combina- 

tion of cn,C66 (more accurate expression for V c will be 
given below, see Eq. (25)). A rough estimate of the 
characteristic velocity V c can be done through the Debye 
frequency, V c ~ ujjj ■ £. The Debye frequency appears 
because of periodicity of the dispersion law and the co- 
herence length indicates that we take into consideration 
all q± wave vectors from \/a v to l/£. (Note V c is much 
smaller that the characteristic phase velocity of collective 
modes, V P h ~ w D ■ a v .) 

The integral in ( p!q ) breaks into the sum over differ- 
ent g. In the framework of a macroscopical approach 
the above derived formula (|l8|) can be presented with co- 
efficients C\ , C2 and A that acquire macroscopic form: 
Ci/al = c n q\ + c^ql, C 2 /a 2 v = c 66 q 2 ± + caq 2 z as were 
defined above and A = (en — c^)q 2 a 2 , C\ — Ci — A = 
(en — CQo)q 2 a 2 . Then, the last bracket in the numerator, 
after replacement q — > q, is equal to zero, that simplifies 
integral calculation. All q vectors those are not included 



in G\ , C2 or A should be replaced by corresponding g vec- 
tors then one can substitute in ( |lq ) both g x — {2n /a v ) n 
instead of q x and g± — (2n/a v ) (n 2 + m 2 ) 1 / 2 instead of 

The sum over m and n may be replaced by the inte- 
gral as we consider large m and n, m,n < a v /£ and we 
integrate over q only a small region near each vector of re- 
ciprocal lattice, \q\ < V/V c , and considering large enough 
values of g till l/£. The expression for the drag force per 
single vortex in the limit of small viscosity, H 3> /3, and 
small velocities, V <C V c , finally is 



F = 1H V 1 ' 2 , (21) 

_ 2na 2 C imp JVH \3/V2, c u ~ c 6 e 

3V c n c 44C 6 6 a^^ 3 / 2 [(cn/c 66 ) 1/4 , en > c 66 , 

where J = J^ 3 ^ 2 dg x dg y f 2 (g±)g x ^ 2 g 2 ^ is the constant of 
the order of unity. 

In the limit of high viscosity, f} 3> H, and V — ► the 
similar estimation demonstrates that again only small 
areas near <f are important and calculations yield 



F = 7/3^ 1/2 ,7/3 = 



\/2TWpCjmp J\f$ 



1 

C66 



1 
Cll 



(22) 



where J is the same constant as in ([2~i]). 

So for small velocities, <C V c , and any relation be- 
tween H and (3 the drag force is F(V) = 7(_ff, /3)v / T 7 , 
where the coefficient j(H, 0) in limiting cases is given by 

©ID- 

When V ~ y c one has to take into account not only 
small regions of the wave vectors near g, but perform 
integration over all \q\ < 1/^ as well. But this is an 
extremely complicated problem and we are not able to 
solve it even numerically as the detailed behaviour of 
Cik at large q is not known. Far above the critical ve- 
locity, V ^> V c , one can significantly simplify the prob- 
lem. In this case, it is possible to neglect contribution 
of the terms like (en, ce^q 2 ^ having the order of the 
value ujo to the denominator in (|l8|), whereas the terms 
C440 2 , /3q x V, Hq x V should be taken into consideration. 
Then, the calculations give at V > V c 



F = 



V 



Jciia^ 1 / 2 



(23) 



where A(H,0) = 1V2H at H > /3 and 
A(H,0) = l/(2y/j3) at H < /?, the constant J' = 
J £, 1 ^ 2 dq x dq v f 2 {q±_)qV 2 q 2 _ ~ J. Below the difference be- 
tween J and J' will be disregarded. The dependence 
F oc \j\/V corresponds to the single vortex limit^l if one 



would replace 044 a 2 



It is not surprising because in 



the limit V ^> V c the only contribution the q z wave vector 
to u>(q) must be taken into account. This is exactly the 
case of the single vortex, for which q x and q y components 
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of the momentum q are absorbed by the vortex, while el- 
ementary excitations propagating along the vortex axis 
acquire the momentum q z and the energy q x V, that does 
not restrict the integration area over q (in contrast to 
© or ©). 

Thus, the velocity dependence of the additional drag 
force, which results from the vortex-defects interaction 
can be described as F = 7 (if, P)W, V < V c , and 
F = rj(H , (3) / y/V , V > Vc, where the coefficients 7 and 
rj are derived in Eqs. ( |2lj - |23j ). To describe the total fric- 
tion force the bare contribution F oare = —f3V should be 
concerned. The total friction force acting on the mov- 
ing vortex lattice having the different velocity depen- 
dence for small and large velocities can be presented as 
F = — (V /V)F(V), with the interpolating equation for 
F(V) 



F(V) 



1 + V/V c 



(24) 



where we use V c = 77 (if, f3)/j(H, (3) as a quantity char- 
acteristic of the problem. The schematic picture of the 
F(V) is presented in Fig. 3. 

This value does not contradict to estimates provided 
above in which only the order of V c has been considered. 
From d2l]-p3"|) one can extract the characteristic velocity 
V c in limiting cases 



V B = 



2V2 



(3/if) (c llC y 1/4 , 
. cnc 66 /[/3(cn + c 66 )], 



if > f3, Cn ~ C 6 6, 
if > /3,Cu > c 66 , 
/?> if. 

(25) 



In the problem one more characteristic velocity 
emerges, Vt rans , at which the bare dissipation starts 
to be dominant compared to the additional one, i.e. 
[3V > jhW/ (1 + V/V c ) at V > Vtrans- We would like 
to emphasize that V tra ns is determined as a velocity be- 
low that the additional drag force prevails over the bare 
drag force, whereas V c is defined only by the additional 
drag force behaviour. In virtue of this, V c and Vtrans 
depend on the problem parameters in different manners. 
As it is seen from ( |2"5| ) , V c does not depend on the lattice- 
defects interaction intensity, i.e. does not comprise Ci mp 
or a. On the contrary, Vtrans is directly correlated with 
the intensity a. Because of the nonmonotonic velocity 
dependence of the additional drag force, Vtrans exists for 
any as small as one likes value of Cj mp a 2 and only the 
character of the total drag force, F(V), is changed. 

If the lattice-defects interaction intensity is small, 
while the bare constant f3 is large, then Vtrans falls on 
the increasing part of F(V) and V tra ns — (7 / ' fif ■ It cor- 
responds to the case Vtrans -C V c and the velocity V c does 
not manifest itself at all: F(V) oc W at V -C Vtrans and 

F(V) CX V at V > Vtrans- _ 

In the opposite case Vt ra ns 2> V c , F(V) cx W at 

V < V c and F{V) OC 1/y/V at V c « V « Vtrans, 



while the linear asymptotic F(V) oc V appears only at 

V > Vtrans > V c and Vtrans = (ry//3) 2/3 . Obviously 
that namely in the latter case the nonmonotonic F(V) 
dependence may be realised. 

The concrete values of this velocity and the relation 
between V c and V tra ns will be estimated below. 



III. FORCED MOTION OF THE LATTICE 

In the previous section we have employed a model of 
a rather weak random force induced by defects acting on 
the lattice. Now we proceed to description of the vortex 
lattice dynamics within a somewhat broader macroscop- 
ical approach. Having use the macroscopical form of the 
average friction force (|2~i| ) instead of the random vortex- 
defects force considered in (^) we are able to reformulate 
the vortex lattice dynamics problem on the macroscopi- 
cal basis. 

Consider the displacement vector of vortices in the 
lattice U = U(r, t) assuming that U changes insignif- 
icantly over the scale of the lattice constant. We use 
the trivial averaging of the Lagrangian (|l|) with the dy- 
namical part presented as in m) and with substitutions 
£;(•••) -> / dxdy/a 2 v {...) and ui(z,t) -> U(f,t). The en- 
ergy of lattice deformation is described by the elasticity 
theory (J3J) . We assume the presence of the external force, 
F e , while the effective drag force is determined by Eq.([2^) 

with the change V — * U. 

The equation of motion takes the similar form as (|^): 



if 



(z x ry) 



ff F(\U\) 



dW{U} 
dU 



= F e 



(26) 



Using the equation (2^)as the base we consider macro- 
scopical dynamics of the vortex lattice, in particular, its 
stability against small perturbations. We present trans- 

lational lattice velocity, U = V, as a sum of the steady 

velocity of the lattice, Vq, and the small correction u, U 
= Vo+ii(r,t). 

First we seek the solution for the lattice steady veloc- 
ity Vq- It allows us revealing F e (V) dependence, which 
is merely current-voltage (I — V) characteristic of the su- 
perconductor, see Refjl^. Because of the Hall term in 
( ptf ) the steady velocity is convenient to express via com- 
ponents Vq = V»e 



V±e±, where eji = F e 



\F e . 



(z x eli). To describe the stationary flow of the vortex 

lattice with velocity Vq a couple of equations of motion 
are found, 



^0 



-HVw 



F(V )Vx 
V 



= 0. (27) 



From these equations one can derive the components 
of stationary velocity V± and Vji via Vq: 
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V\ = 



F e /H 



1 + (F/HVo) 



FF e /H 2 V 
1 + (F/HVc 



(28) 



For simplicity we use F' 
F(V Q ). With the help of | 
given as 



[dF(V)/dV]\ v=Vg , F = 
the Hall angle could be 



tana = 



V± 
V» 



HVp 

F(v y 



(29) 



So if the drag force is a linear function of the transla- 
tional velocity Vq, F = f3Vo (the usual case), then tana 
does not depend on the velocity and the Hall angle tan- 
gent is merely the ratio of the Hall constant and the 
dissipation constant j3. However, by considering the spe- 
cific drag coefficient in Eq.(^J) a striking feature appears, 
namely a velocity dependence of the Hall angle. In the 
limit V — > tana = \/Vo/j(f3,H), and the higher the 
lattice velocity the Hall angle larger, then it is saturated 
above the characteristic velocity V trans . 

By solving the E q.(|2q ) with the taking into account the 
drag force from Eq. (|24|) the dependence of F e (V) acquires 
the following form: 



V?H 2 



1 + Vo/Vc 



= Fi 



(30) 



The Eq.(|3^) has a simple form of the vector sum of two 
forces: the Hall force (the first term on the left hand side) 
and the drag force (the second term). The differential 
mobility related to the external force is 



dVo_ 
dF R 



V H 2 + F{Vv)F'(V y 



(31) 



If one puts H = in (|3l| ) then the differential mobility 
dV /dF e coincides with (dF/dVo)- 1 . But for H ^ 0, 
especially at H ^> F(V)/V these two quantities differ 
substantially. 

The negative value of the differential mobility is usu- 
ally considered as a condition (both necessary and suf- 
ficient) of the instability appearance. Nevertheless, it is 
quite possible that dVo/dF < 0, whereas dVo/dF e > 0, 
see @. Therefore, one has to reconsider the instabil- 
ity criterion for the case of macroscopic dynamics in the 
presence of gyroforce. This problem is of interest not 
only Abrikosov lattice dynamics, but it is also important 
for dynamics of any system having gyroforce and a non- 
monotonic velocity dependence of the drag force, F(V) 
(for instance, vortices in He 4 , Bloch lines, etc.) 

To investigate the stability of the vortex lattice we 
take U — V$t + u(r, t) and find the solution in the form 
u(r,t) oc w 9 exp(Ai + iqr), where A defines the charac- 
ter of time evolution of small deviations from the steady 
motion. If the real part of A is negative, the value of 
— ReA > serves as a damping coefficient. Otherwise, if 



the real part of A is positive for some values of q, cor- 
responding small deviations cause the motion instability 
with the instability increment ReA > 0. 

For the small correction u the couple of equations of 
motion arc 



u q . x \AF'{V) + C xx ] + u q , y (C xy - AH) = 0; (32) 
u q , x (C xy + AH) + u q . y [AF{V)/V + C yy ] = 0, 

where Cik = Cik(q) are the coefficients used in (|l5|), 
strictly speaking, their presentation in the longwave 
limit, V is the translational velocity of the lattice (for 
simplicity we omit here and futher the index 0) . It is more 
convenient to use in Eq.(|32|) e x , e y instead of en, e±. It 
is necessary to stress that this equation contains F(V) 
and F'(V) = dF(V)/dV, where F is the total drag force, 
rather than differential mobility dF e /dV . As it will be 
demonstrated below, this is a very essential point which 
leads to the fact that the instability condition is not re- 
lated directly to the usual one, dF e /dV < 0. 

The characteristic equation for A takes the form: 



A 2 [F'(F/V) + H 2 ] +uj 2 (q)H 2 
+ A [F'Cyy + [F/V)C XX ] = 0, 



(33) 



where we used the equation C xx C yy — C xy = uj 2 (q)H 2 , 
u>(q) is the dispersion law for free vortex oscillations (||) 
obtained above. The equation (p3|) possess two solutions: 



A = -r(g)±*n(«), 



where 



r(«) 
n 2 ( q ) 



(F/V)C X 



C, F' 

^j yy r 



H 2 + (F/V)F' 
H 2 L0 2 (q) 



(34) 



(35) 



H 2 + (F/V)F' 



In order to get the instability condition we analyse 
Eqs.(||, HI)- If F' > 0, then r > and small oscillations 
decay for all values of q. These are small damped oscilla- 
tions with the frequency u) w £l(q) at il <C T, otherwise, 
for £1 ^> r, this is the case of overdamped dynamics of 
u. Thus, the inequality F' < is the necessary condition 
for the lattice instability. 

But if F' < 0, it is not sufficient to induce the lattice 
instability. The further analysis of Eqs.(|34|, |35| ) reveals 
two different and independent sufficient conditions im- 
posed by the instability requirement ReA > 0. 

The first condition is H 2 + (F/V)F' < or 



-F'(F/V) > H 2 



(36) 



This requirement is equivalent to the inequality 
dF e /dV < 0, i.e. the rise of the negative differential 
mobility, in a literal sense. When this condition is sat- 
isfied the values of A in (|34|) are real and the instability 
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appears at any q ^ 0. [If one considers the inertial term 
in the equation of motion (ffif), then even the uniform 
deviations U(t) from Uq (dynamics with q = 0) become 
unstable when dF e /dV < 0.] 

Hence, the fulfillment of the condition dF e /dV < ap- 
parently indicates the rise of the instability of the lattice 
translational motion. But this condition is substantially 
more strict than F'(V) < (see Fig. 3), and as we will 
demonstrate in the next Section never will be satisfied in 
our model (this condition, though, could be met in other 
models with dF/dV < 0). 

The second condition corresponds to the inequality: 



-F'Cyy > (F/V)C XX 



(37) 



The condition ( p7\ ) is one more sufficient condition of 
the instability independent of (p6|). Depending on the 
strength of requirements ( ^6p^ ) one of them should be 
fulfilled earlier than the other. We will demonstrate be- 
low that the condition ( |37| ) will be met in real supercon- 
ductors unlike the requirement (|3^). 



IV. DISCUSSION 

In the previous section we have established a set of con- 
ditions of the vortex lattice instability which comprises 
both the necessary and sufficient requirements. Let us es- 
timate the least strict condition, the necessary condition 
F'(V) < 0. By differentiating one obtains 



7 (V/V c - 1) 



2vv (i + v/v c y 



>ft. 



(38) 



The left hand side of the inequality (|38|) as a function 
of the variable V/V c has the maximal value O.OAj/vVc 
at V/V c ~ 2. Using V c = rj/j one can rewrite (|3§) as 



ft < 0.04 • 7 3/2 /\A? 



(39) 



The parameters r\ and 7 can be taken from the limit- 
ing formulae ( pjlf]23] ) and expressed through microscopic 
characteristics of the superconductor (see, e.g., Refs. [|, 
|l6| ): the elastic modulii Cu = C44 = B 2 /4tt, cqq — 
B$o/(8ir\) 2 , where B is the magnetic induction, $0 
is the flux quanta, A is the penetration depth. The 
vortex lattice constant is determined by magnetic in- 
duction a 2 = (2/V3) ($o/B). The Hall constant is 
H = &oen/c, where e is the electron charge, c is the speed 
of light, n is the density of superconducting electrons, 
n = mc 2 1 '4ir\ 2 e 2 with tliCL-plectron mass, m. The param- 
eter a can be written as! a = a 3 (H 2 m /8Tr) (AT C /T C ), 
where H cm is the superconducting thermodynamic criti- 
cal field, AT C /T C is a relative suppression of the critical 
temperature nearly a defect, a 3 is the volume of the de- 
fect (a is of the order of the interatomic distance). As 
a result, in the case of the small dissipation, the Eq.(|39|) 



takes the following form: 

ft ^ 11 s\ f^T c \ 2 a 3 \ 3 [H cl \ p , s 

where H c i = $0/ 47rA 2 is the quantity of the order of the 
lower critical field and p depends on (ch/cqq) 



19/8, en > c 66 , 
11/4, cu ~ C 6 6. 



(41) 



If one would take the numerical values £>f A and £, as 
for typical high T c superconductor YBCO:El A = 10~ 5 cm, 



£ = 10~ 7 cm, a = 5 • 10~ 8 cm, AT C /T C 
1Q- 3 then 



ft 



l,C, 



imp 1 - 



{H cl /Bf > J£ 



(42) 



This condition is met at magnetic fields close to the 
lower magnetic field, H c i, and the microscopical param- 
eters used above for YBCO compound. The maximal 
values of (3 = ^ 2 j2n^ 2 c 2 p n , where p n is the normal con- 
ductivity of the superconductor. So the maximal value 
of the ratio ft/H is ft/H = (1/ p n )(\/£) 2 (h/27rmc 2 ) and 
at p n = 5 • 10~ 16 s (the typical value for YBCO) it is 
equal approximately to 10~ 2 . Note this is just the upper 
limit of ft/H and for real superconductors this value can 
be even smaller. Thus, for the case of low viscosity the 
instability necessary condition may be fulfilled at high 
enough values of magnetic fields. 

In the case of high viscosity, H <§C ft, the necessary 
condition can be satisfied with difficulty as lift) contains 
ft and the right hand side of the inequality ( pq ) becomes 
unity (in conrast to the case H ^> ft where the small 
parameter ft/H takes place). 

Now we study the first sufficient condition 
—F'(F/V) > H 2 . We carry out the calculation in 
the similar manner as for (B9) and obtain (neglecting 
the bare dissipation as H » (3) that 0.1 • "f 3 ^ 2 /^/fj > H. 
Apparently this condition cannot be satisfied since 
though the same inequality as ( p2[ ) is valid for this 
condition, but on the right hand side instead the small 
parameter ft/H unity is found. Then, we consider the 
second sufficient condition —F'j (F/V) > C xx /C yy 
which acquires the following form: 



7 ((C yy ~ 2C XX ) V/V C - (Cyy + 2C XX )) 
(Cyy + C xx ) (1 + V/V c f 



2VV 



>ft. (43) 



It is seen that the inequality (|43|) has an obvious re- 
semblance with Eq. (|38|) . The difference is caused by q- 



dependent coefficients containing C. 
For 2C X . 

cu — cqq, the left hand side of (|43| 



and Cy y in Eq. (k 



> Cy y , which is realised, for instance, when 



is negative. There- 
fore, for some values of C xx , C yy the condition ( f43|) can- 
not be fulfilled and no instability appears. However, by 
using the fact that c 66 — (cn/4)(iJ cl /_B) <g; c n and 
choosing an appropriate g* it can be satisfied. The simple 
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analysis shows that when cu > cqq the most favorable 
case is the wave vectors along the y-axis, q\\e y . Then, the 
inequality (j43|) takes the following form: 



7 ((ci i-2c 66 )V/V c -(cii+2c fe6 )) 
(en + c 66 ) (1 + V/V c f 



2VV 



> 0- (44) 



The maximal value of left hand side of ( p4[ ) is attained 
at 3 (cu — 2c 6e ) (V/V c ) = (3cn + 4c 6 6 + 2c), where c = 
a/Sc^ + 6c n c 66 + c| 6 . It is equal to 



7 3 / 2 (3V3/4) (en - 2c 66 ) 5/2 (c - c 66 ) 

(cn + c 66 ) V3cn + 4c 66 + 2c (3cix - c 66 + c) 



Thus, the sufficient condition 
necessary one given by 
cu and C66- In particu 



(45) 

is more strict than 
owing to the multipliers with 
ar, it can be met only at cu > 
2c66- However, in the case of interest cu > 10 cqq, the 
necessary and sufficient conditions are close each other 
and the sufficient requirement (|44| ) can be satisfied as 
well. As the instability is expected to arise when H ^ 
(3 and V » 2V C , one can estimate V c for this case. In 
accordance with Eq.(^5|) and with the use of the above 
presented microscopical parameters for YBCO V c can be 
presented as: 



V c « 10" 2 • (B[Oe]) 3/2 , 



(46) 



and at B = 10 3 Oe, V c w 300 cm/s. Thus, the instability 
may appear under the rather strict condition ( |42"| ) and at 
the velocities of the order of 2V C ~ 10 3 sm/s. Hence, both 
the necessary and sufficient conditions can be fulfilled at 
low viscosity, (3 <C H, relatively small magnetic fields and 
quite high velocities of the vortex lattice motion. 

At small velocities another effect, namely a nonlinear 
I- V characteristic of the superconductor emerges. This 
effect does not impose limits for superconductor parame- 
ters. It can be observed at (3 > H and at strong magnetic 
fileds, when V c 3> Vtrans and no instability arises. 

At V c ^S> Vtrans, F(V) dependence is monotonic and 
contains only part with F cx s/V, V < Vtrans, F cx 
V, V > V^ans ■ In this case, the quadratic I- V character- 
istic is realised at V < Vtrans ■ This effect can be found at 
very small defect concentration and rather dense lattice. 
For instance, if Ci mp a 3 — 10~ 6 , Vtrans ~ 10 cm/s when 
H > (3 and Vtrans « 70 cm/s for H < (3 at B = 10 3 Oe. 



V. CONCLUSION 

We have demonstrated that due to the interaction of 
weak defects with the Abrikosov vortex lattice the addi- 
tional channel of dissipation appears, which contributes 
to the total dissipation of the moving lattice. As a re- 
sult the additional friction force has a non-linear depen- 
dence of the translational lattice velocity, F cx \fV for 
both weak bare dissipation, H ^> f3, and high viscosity 
H (3. This channel of dissipation turns out to be an es- 
sential source of the system dissipation and even prevails 
over the bare dissipation, (3V, below the characteristic 

velocity Vtrans- 

The nonlinear character of the additional drag force 
results in a nonlinear current-voltage characteristic of 
the superconductor. Besides, the tangent of the Hall an- 
gle tana appears to be the translational lattice velocity 
V) dependent up to the characteristic velocity Vtrans, h, 
above which it is saturated. 

Studying the stability of the vortex lattice against uni- 
form and nonuniform perturbations we have found sev- 
eral salient features. 

The instability of the vortex lattice motion is found to 
occur in the limit of small viscosity, H ^> (3. We have 
revealed that the negative differential mobility, dVo/dF e , 
is not apriori the lattice motion instability requirementl. 
Due to the Hall component in the lattice motion, the 
value [dF(V)/dV]~ 1 does not coincide with the differen- 
tial mobility dV/dF e . The estimates suggest that condi- 
tion dV/dF e < can be met with difficulty in contrast 
to dV/dF < 0. The instability of the lattice motion re- 
garding nonuniform perturbations has been found when 
the value of (-dF/dV) is greater than the combination of 
the elastic modulii cu and c^e, while the more strict re- 
quirement dV/dF e < is still not fulfilled. It is worth to 
note that the vortex lattice instability concerned in this 
paper has different origin than instabilities predicted in 
Ref. 0. 
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FIG. 1: The sketch of the vortex displacement ui(z,t). 




FIG. 2: The schematic representation of the dispersion law 
of the vortex lattice. The dotted line indicates the Debye 
frequency, loo- 
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FIG. 3: The velocity dependencies of the drag force F(V) 
at B — 3-ffci for different values of the ratio (5/H (its value 
is shown near each curve). The dotted line parallel to the 
x-axis marks the point on the curve with the critical ratio 
13/ H = 0.0074 where the nonmonotonic character of F(V) 
starts to develop. The dashed line describes the relation of 
the external force, F e and the velocity V , at very small value 
of the ratio /3/H — 0.002. As was noted in the text the 
behaviour of this curve is very far from nonmonotonic and 
can be hardly distinguished from the y-axis. 



